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ABSTRACT 

The discovery of multiple evolutionary sequences has challenged the paradigm that globular clusters 
(GCs) host simple stellar populations. In addition, spectroscopic studies of GCs show a spread in light- 
, element abundances, suggesting that multiple sequences can be formed from gaseous ejecta processed 

' in evolved cluster stars. If multiple sequences originate from within GCs, then it should be determined 

how such stellar systems retain gas, form new stars within them and subsequently evolve. Here we 
expand upon previous studies and carry out hydrodynamical simulations that explore a wide range of 
cluster masses, compactness, metallicities and stellar age combinations in order to determine the ideal 
conditions for gas retention. We find that up to 6.4% of the mass of the star cluster can be made up 
of retained stellar wind gas at the time star formation is triggered. However, we show that multiple 
(N . episodes of star formation can take place during the lifetime of a star cluster in particular for times 

> 1 Gyr, thus leading to a sizable enhancement in the total number of new stars. The fact that this 
favorable star formation time interval coincides with the asymptotic giant branch (ACB) phase seems 
to give further credence to the idea that, at least in some GCs, there are stars which have formed 
from material processed by a previous generation of stars. The ability of extended heating sources, 
such as pulsar outflows or accretion onto compact objects, to hamper gas retention is illustrated via 
a simple numerical treatment. 

Subject headings: globular clusters:general, hydrodynamics, accretion 



(N 



e3 

^3 



O 



O 

o 



>< 



O 

^ ' 1. INTRODUCTION 

The idea that stars in globular clusters (GCs) are a coeval population has been a long disputed topic wi th significant 
discussion devoted to the presence (or lack thereof) of multiple main sequences and subgiant branches (see lPiottoll2009L 
for a review). However, recent observational evidence has demonstrated that these multiple stellar pop ulations (MSPs) 
^ ' are not only ubiquitous but a lso make up a sizable fraction (~40-60%) of the total stellar mass (Cal oi fc D'Antonal 
[ 120071 : iD'Antona fc Caloill2008l). Spectrosc opic observations of MSPs in GCs have also unveiled differences in their light 
element abundances (see e.g. iGratton et al. 2004: de Silva ct al. 2009; M artell fc Sm ith 2009). These anomalies can 
only be produced if these MSPs were formed from material processed at temperatures T > lO^K (jD'Ercole et al.l[2008l : 
IRenzini fc Voiil[T98lHVmtura fc D'Antonal I2008bn . 
^^O . To explain the presence of MSPs with the observed abundance anomalies the followin g timeline has been identified 
O • ()Conrov fc Spergel|[20Tll: iCottrell fc Da Costal fT98lt iSmithI [19871: iCarretta et aIll2010bD : after the first generation of 
[ stars is formed and intracluster gas is expelled by supernovae, matter processed in the interiors of AGB stars is 
returned to the ISM through winds and after several 100 Myrs a second generation of stars is formed from a mixture 
of stellar-processed material and captured ISM gas. Some of the issues at the forefront of attention include the type 
and age of stars able to process gas at the required internal temperatures while at the same time efficiently returning 
material into GCs to form a sizable mass fraction of new stars. We address both of these issues here. 

The effective retention of stella r wind inateri al to create MSPs of comparable mass presents a challenge to the current 
understanding of GC formation (!Conrovll2011[ ). Accounting for the inflow and mixing of pristine ISM material slightly 
lowers the retained mass requirements, however the original GCs still need to be significantly more massive than those 
curre ntly observed (jNaiman et all 120111 : iPflamm-Altenburg fc Kroupal [2009t iD'Ercole et"all 120101 : iConrov fc Spergell 
1201 It) . D espite their relevance, stellar wind gas retention has so far only b een studied for a very restricted number of 
systems ()D'Ercole et al.|[200l [20Tol: lConrovll201l"l: IConrov fc SpCTgei|[20TTI ). Here we expand upon these studies with 



hydrodynamical simulations which explore the possible combinations of cluster mass, compactness and stellar age to 
determine the optimum parameter space for maximal gas retention. 

2. NUMERICAL METHODS AND INITIAL SETUP 

2.1. Hydrodynamics 

To examine the ability of clusters to effectively retain the winds emanating from their evolving stellar members, we 
simulate mass and energy injection in isolated core potentials under the assumption of spherical symmetry (Ouataer^ 
[200l iHuevotl-Zahuantitla et al.ll2010l ). The one-dimen sional hydrodynam ical equations are solved using FLASH, a 
parallel, adaptive mesh refinement hydrodynamical code (jFrvxell et al.ll2000l ). The winds from the closely packed stellar 
members are assumed to shock and thermalize such that density and energy contributions can be treated as source 
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term s in the hydrodynamical equation s. In spherical symmetry, these equations can be written as (|Holzer fc AxfordI 
flfza iHuevotl-'Zahuantitla et al.ll2010l) : 

f + ;^|:(^-^)-»»W (1) 

du du IdP ^ ^ 

l + ^^M + ^l^^'/eM-QM (3) 

where P — P{r), p = p{r), u = u{r) and e = e(r) are the gas density, pressure, radial velocity and internal energy 
density, respectively. Here, Q{r) = ni{r)ne{r)A{T, Z) is the cooling rate for gas with ion and electron number densities 
ni{r) and ne{r), and A(T, Z) is the cooling function for gas at temperature T and metallicity Z. We use A(r, Z) from 
IGnat fc Ster^AeFi (|20^ for T > 10^ K and Da karno fc McCrav (1972) for 10 < T < lO"* K. 

The qm{r) and ge(r) terms in equations (P)-© are used here to mediate the total rate of mass and energy injection. 
For N stars each with average mass loss rate (M) and wind energy injection rate \{M){v1,) we have M^, total = 
N{M) = J 47rr2g„(r)dr and i^w.totai = ^N{M){vl) = J 47rr^q^{r)dr, such that q^{r) = ig^(r)(?;2 ). For simplicity, 
we neglect the effects of mass segregation and assume qm(r) oc n{r) cx r^"^-^ {f^'lT^- ^^^'^ neglect the effects 

of external mass accumula tion, which can be an important source of gas when clusters reside in cool, dense ISM gas 
aiman eral|[2009l . [20Tl[) . While specific heatin g sources such as photoionization or supernova are not explicitly 

included, the overall effects of additional energy injection in the cluster are discussed in jjH 

Th e stellar cluster gravitational potentials are modeled as Plummcr models (jBriins et al. l2009H Pfl amm-Altenburg fc Kroupal 

[20091) 

e^^^ ,4) 



for a given total mass, Mc, and velocity dispersion, = (3''/''/ V2) GMcjr^. In addition to following the dynamics 
of the gas under the influence of $g, the self gravity of the gas is computed using FLASH'S multipole gravity module. 

In cases where catastrophic cooling occurs, a sub-gri d model for star formation is implemented to prevent premature 
collapse. Our cooling prescription is modified following iTruelove et al.l (|1997| ) by turning off energy losses and allowing 
the gas to evolve adiabatically when the Jean's length of the collapsing gas structure is comparable (or smaller than) 
to the cell size of our simulation. Such a method approximates the transition of an isothermally collapsing cloud 
to an optically thick, adiabatically evolving protostellar cluster while forgoing the computationally expensive three 
dimensional radiative transfer calculations required to treat this problem accurately. Because we do not include an 
explicit star formation prescription, we allow such unstable regions to evolve for a few sound crossing times before 
halting the simulation. 

2.2. Stellar Evolution 

To est imate the net mass lo s s and mean thermal velocities of the colliding winds, we first follow the formalism devel- 
oped bv lPoolev fc RappaportI ()2006[ ) where the integrated wind kinetic energy and mass loss of the stellar population's 
turn off star is used as a proxy for the average wind velocity and mass loss rate: 

(«w,to> « (5) 

Here AEk = ^ J^^ Mtov'^ to'^t and AM = J*^ Mtodt are the kinetic energy and mass loss input rates integrated over 
the lifetime of the turn-off stars, At = t2 — ti. 

The effects of the additional input of mass and kinetic energy by stars with M^, < M^-^to, which are neglected in the 
formalism we just outlined, can be included by convolving the above definitions of th e average mas s loss r ate and wind 
velocity with an initial mass function (IMF) and a star formation history. Following iKroupa et al.l ()2011| ). the average 
number of stars in a mass interval [M^,, M^, + dM^] evolving between [t, t+dt] is given b y dN = ({Mi,,t)b{t)dMi,dt where 
(^{Mi,,t) is the IMF, which we assume is accurately described by the iKroupal ()2001[ ) IMF, and b{t) is the normalized 
star formation history. For a non-evolving IMF, C(M*,t) = ({Mi,), with a mass distribution extending from Mj,_i to 

Mi^^H, the normalized star formation history can be written as jj='s'='^^*'^' b{t)dt = 1, where iagc(-^^*) is the lifetime 
of a star of a given M^,. Examples of b{t) include a constant star formation rate, b{t) = 1, or a coeval population, 
b{t) = 5{t — to), with all stars formed at to- 
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The average mass {AM(ti)) and kinetic energy {AEk{U)) injected into the cluster environment at a time ti, by a 
population of stars with Af^^e [M^,^L, Af^^jj] and whose birth rate is regulated by b{t), are then given by 

{AM{ti))= b{t) (:{M,)M{M^,t)dA'hdt (7) 

J to Jm,,i, 



and 



1 



{AEk{U))^- h{t) CiM.)M{M,,t)vi{M^,t)dM,dt, (8) 

^ J to "'m,,l 

respectively. Stars that have lifetimes tago(Af*) < ti forsake the stellar population unless they were born during the 
most recent star formation time interval [U — ta.ge{Mi,),ti]. For this reason, to is set to max[0,ti — tugdM-i,)]. This 
formalism allows for equations ([5]) and ^ to be cast into a more general form 



(AMiU)) 1 



and 



(«^(iO) = ^||§^ = (Af(t.)>"'y^ bit)J^^ aM.)M{A'U,t)vliM.,t)dM.dt. (10) 

In what follows, we assume the stellar population initially residing in star clusters is coeval such that b{t) = (5(0), 
which provides an accurate description for clusters with ti > 10 Myrs. 

To co mpute the individua l mass loss rates M{Mi,, t) and wind velocities v^iMi,, t) we use the MESA stellar evolution 
models (jPaxton et al.|[20TTl ). We use this code to follow the evolution of a grid of stellar models with Mi, — 0.1 Mq 
- 20 Mq from ZAMS to either the white dwarf stage, end of the AGB phase, or compact ob ject creat i on. T he mass 
loss rates during the crucial AGB phase are estimated here usi ng the wind prescrip tion of iBloeckeil (|1995l ) with a 
normalization of rjB — 0.04 consistent with LMC measurements ([Ventura et al.lf2000( ). We fix Z = I/IOZ© for most 
models although the effects of changing the metallicity are discussed in 21 

3. STELLAR WIND RETENTION IN STAR CLUSTERS 

The mass injection properties, characterized here by the mass loss rates and wind velocities emanating from a coeval 
population of stars, change dramatically as the population evolves. Figure [T] shows the evolution of the average stellar 
wind parameters. As the wind velocities and stellar mass loss rates change, so does the ability of a cluster potential 
to retain the shocked stellar wind gas. Note that the differences between the turn off mass and population averaged 
prescriptions manifests itself predominately in the values of the wind velocity. 

A young star cluster, corresponding to time A in Figure [U with its high average wind velocity cannot retain gas 
effectively if <C Vw This is clearly seen in Figure [2] which shows the properties of the shocked stellar wind gas 
confined to a cluster potential characterized by ay = 50 km/s and Mc = 10^ Mq for both turn off mass and population 
averaged prescriptions {red curves). The fiow in this case was evolved for ti ~ tA- A young cluster is thus only likely to 
retain a small quantity of high temperature gas in its inner region, while further away gas is blown out of the system in 
a wind. As the cluster members evolve, a dramatic decrease in the average wind velocity occurs due to the dominant 
contribution of the slow, dense AGB winds to the overall mass injection, corresponding to time B in Figure [T] For the 
star cluster, here assumed to be characterized by a non-evolving gravitational potential, the AGB contribution results 
in stellar wind injection parameters that favor significant mass retention, as shown in Figure [T] (6^acA; curves). As the 
stellar wind gas shocks and cools, mass is efficiently accumulated until the central region becomes Jeans unstable, thus 
triggering star formation. In this case, the gas flow within the cluster is evolved until t — tgf, defined here as the sound 
crossing time within the (adequately resolved) Jeans unstable region. As the cluster ages, remaining stars no longer 
pass through an extended thermally pulsing AGB phase. This results in a dramatic increase in average wind velocity 
and decrease in average mass loss rates as seen at time C. As a result, the cluster potential is unable to effectively 
retain the emanating gas and the wind material ffows out of the cluster almost unrestrained. 

The ability for a star cluster to retain wind ejecta depends not only on the evolutionary stage of its members but 
also on their spatial distribution. To illustrate this, in Figure [3] we show the properties of the shocked stellar wind 
gas confined to a shallow potential characterized by ay = 14 km/s and Mc = 4.5 x 10^ Af ^; ), parameters though t to 
accurately represent the current stellar mass distribution in Ml 5 (jMcNamara et all I2Q04I : iGerssen et al.l [20021) . A 
shallow gravitational potential is thus unable to retain a significant amount of stellar wind materi al, even during the 
AGB phase. The observational constraint on the allowed gas density in M15 (jKnapp et al.l[T996[ ) is consistent with 
our predictions for the present state of the gas in this system. 

The dramatic difference observed in Figures [2] and [3] between the stellar wind mass accumulated in cluster potentials 
of varying properties motivates our study to compare results obtained with different velocity dispersions and cluster 
masses. To facilitate comparisons, we first systematically vary the cluster velocity dispersion for a fixed total mass 
Mc = 10^ Mq, assuming the turn off mass wind velocity and mass loss rate are representative of the entire population. 
The left panel of Figure |4] gives the amount of accumulated stellar wind mass in units of Mc for a range of cluster 
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velocity dispersions and different evolutionary stages (t^) of the stellar members. By looking at the shaded regions in 
Figure 21 the reader can identify the cluster velocity dispersion and mass combinations for which favorable conditions 
for star formation are satisfied at a given cluster age. 

Cluster potentials with cr„ < 20km/s are not effective at retaining gas in their cores, and gas is blown out of the 
cluster at all times. As the velocity dispersion increases above 30km/s, the cluster core is able to retain the shocked 
and efficiently cooled stellar wind gas, making the central region Jeans unstable before (i.e. tgf <^ U) a significant 
amount of mass is accumulated into the cluster. The largest fraction of mass retained takes place for star clusters with 
(jy « 25 km/s. In this case, the potential is steep enough to retain the shocked and subsequently cooled wind material, 
but shallow enough to force the cold gas to collapse only until a significant amount of mass has been accumulated. 
The contours in the left panel of Figure |4] show that potentials with (t^ > 30 km/s can endure multiple episodes of 
star formation albeit involving less retained gas, while potentials with 25 km/s ^ cr^ < 30 km/s can have only short 
but more intense star formation periods. The largest mass accumulation for a single star formation episode occurs for 
ay « 27 km/s between the ages of 1 — 10 Gyrs, with a mass retention fraction of M ^^cc/Mc « 6.4%. This value is lower 
than the salO% estimated by previous studies (|Conrovll201lHD'Ercole et al.ll2008D . owing to the fact that we used the 
time averaged values of the wind velocity which results in about five times larger kinetic energy injection rates than 
using the instantaneous values. 

The right panel of Figure |4l in conjunction with the left panel, allows us to estimate the amount of stellar wind 
material retained by a star cluster of a given age, mass and velocity dispersion. The conclusions drawn are based 
upon the assumption that the cluster potential properties, whose evolution causes are poorly known, remain relatively 
unaltered for mm[ti,ts{]- The right panel of Figure |4] is self-explanatory - heavier star clusters at a fixed velocity 
dispersion retain more gas. The mass retention fraction is observed to increase slightly with augmenting cluster mass 
as a result of less efficient cooling, which in turn delays the triggering of star formation as mass continues to accumulate. 

Using the average wind velocity and mass loss rate of the turn off mass star to represent the properties of the entire 
population underestimates the kinetic energy input arising from the more numerous, lighter stars (Figure[T]). This gives 
an optimistic value of the wind retention and star formation efficiencies. The effects of using the total kinetic energy 
input from the star cluster to calculate effective mass retention are shown in Figure \5\ Direct comparison with Figure 
|4] shows that although the overall mass accumulation is rather similar, in the population-averaged prescription star 
formation is triggered over a narrower range of velocity dispersions. Not only are higher velocity dispersions required 
to produce efficient star forming models, but the additional kinetic energy inhibits star formation in the largest mass 
accumulation regions. This is because the larger wind velocities keep the gas at higher temperatures, thus quenching 
star formation. This mass is subsequently blown out of the cluster at later times when the mass loss rates decrease and 
the accompanying wind velocities increase. These results clearly indicate that care must be taken when attempting 
to estimate the amount of new stars that can be subsequently formed in star clusters as this strongly depends on the 
assumed mass loss prescriptions, which are uncertain. 

While the assumed wind prescription can alter the total amount of gas retained which can create new stars, the 
contours in Figures |4] and [5] suggest that a stellar cluster with favorable potential parameters, [Mc, (7„], can trigger star 
formation over a large span of time. In reality, an episode of star formation will trigger subsequent supernova outbursts 
which can drive gas away from the central regions of the cluster, deterring further star formation. We simulate the 
effects of supernova explosions on gas retention by adding thermal energy to the central region of a typical cluster 
once star formation has been triggered. To maximize the effects of feedback in our cluster, we assume a large star 
formation efficiency of 90%, creating M^, w 10^ Mq in se cond generation stars for a model from the left panel of Figure 
IHwith a 6.4% retention rate. Following the method of iD'Ercole et al.l (|2008[ ). we add in Esn = lO^^ergss"^ in the 
central regions of our simulation at a rate of 2 x 10~^yr~^. This ha s been scaled from their supernova type la rate of 
2 X 10~^yr~^ for a Mc = 10^ population (jMarcolini et al.l[2006( ). After a period of 100 Myrs the supernova energy 
injection is halted and the gas is once again allowed accumulate in the central regions of the cluster, starting the 
cycle over. Figure [S] shows this cycle of mass accumulation, star formation and supernova feedback for a simulation at 
tage ~ 2000 Myrs for a model with Mc — 10^ Af©, (t^ — 26 km/s, and Macc/A/c = 6.4%. Here, it takes approximately 
2000 Myrs for the cluster to re-trigger star formation once gas has been removed by 100 Myrs of type la supernova 
feedback. 

4. DISCUSSION 

Throughout this work, we have used simplified models to determine how effective a star cluster of a particular age 
is at retaining gas emanating from its stellar members, under the assumption that its potential remains unaltered 
during the simulated phase. We have further assumed a single metallicity for all clusters and have disregarded any 
heating sources besides supernova feedback and the stellar winds themselves. In this section we relax both of these 
assumptions. 

4.1. Metallicity 

Clu ster to cluster variations in light element abundances are commonly observed (j Caldwell et al .11201 It iBeaslev et al.l 
120051 ) . These variations may cause changes in the mass loss histories of the individual stellar members and the cooling 
properties of the shocked gas. As the stellar mass loss prescription are mostly independe nt of metallicit y during the 
evolutionary time periods that are conducive to star formation (0.1 Gyrs < < 100 Gyrs) (lBloeckeiill995f) . we account 
for the effects of varying metallicity solely in the cooling function. Figure!?] shows the effects of changing metallicity for 
a typical star forming cluster, described here by Mc = 10 Mq, (7„ = 50 km/s, and ti = 212 Myrs. As the metallicity is 
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decreased, the cooling becomes less efficient and more material is allowed to flow into the center of the cluster before 
catastrophic cooling occurs. This enables relatively more massive star forming episodes to be triggered. Interestingly, 
when metallicity is decreased beyond Z < 10~^, the cooling becomes weak enough to prevent catastrophic cooling at 
times < ti. These results suggest that the range of cluster parameters over which large central densities will persist 
before catastrophic cooling takes place (Figure^) will depend on the metallicity of the emanating stellar winds, though, 
as illustrated in Figure [7l the differences are not marked. 

4.2. Intercluster Heating Sources 

In addition to altering the cooling curves, the inclusion of additional cluster heating sources may prevent effective 
gas retention in our simulations. We address this problem here by artificially increasing the energy input rate: ge.now — 
(1 + H)q^ = (1 + _ff)J(7m(r)w^. Under this assumption, the additional heating sources follow the potential's stellar 
distribution. Figure [8f shows the effects of the additional heat input in one of our otherwise star forming simulations. 
For H < 2.0, the gas in the simulation still collapses, triggering star formation. For larger values, on the other hand, 
the cluster is unable to effectively retain the gas and, as a result, star formation never ensues. By integrating Qe,ncw 
over the cluster's core for H ~ 2.0, we derive the total energy input rate required to overturn the central mass build 
up, which for this simulation is about 10^^ ergs~^. 

In many cases, the additional energy injection sources might not follow the stellar distribution. As an example, 
let's compare the heat distribution expected from accreting neutron stars under the assumption that the accretion 
feedback is proportional to the Bondi accretion rate: g^^ns oc M c>c p(r)T(r)~^/^ (.Bondi fc H ovle 1944j). Using the 
volume-averaged density and temperature in the cluster core, p « 10^^^ g cm~^ and T « lO'^K, we derive the average 
luminosity of a single, accreting neutron star: Lj^s ~ 10'^^ ergs s^^. This implies that > 100 accreting neutron stars 
are required to reside in the cluster's core in order to significantly offset its cooling properties. However, to accurately 
test this phenomena a multi-dimensional approach would be required as feedback would not necessarily act as a simple 
heating prescription. 

In this work, we have also tried to minimize the effect of external mass inflow by considering star clusters in 
isolation. This is certainly not the case for clusters moving through cold, dense envir onments, as they can potentially 
amass a significant amount of gas from their surroundings ()Naiman et al.|[2C)09ll201l"l ). The stellar winds and exterior 
inflows in s uch clusters could combine to create even larger central density enhancements (j Naiman eFaDiiooi unn 
iConrov fc~ Spergel 2011; Pricstlcv ct al. 2011; Pflamm-Altcnburg & Kroupa 2009). A self consistent treatment of both 
interior and exterior gas accumulation in multi-dimensional simulations, which includes the effects of compact object 
accretion, tidal stripping, photoionization, and pulsar heating will be presented elsewhere.^ 

While previous work has estimated the ability of star clusters to retain stellar winds ([D'Ercole et al.l 120081 120101 : 
IConrovll2011t IConrov fc Spergelll2011h . calculations have so far been restricted to a small range of cluster properties 
and stellar ages. Motivated by this, we have calculated gas retention in star clusters of various ages, stellar mass and 
compactness. In agreement with previous studies we find that before star formation is triggered about <10% of the 
total cluster mass is comprised of retai ned stellar wind gas, naively implyin g that the original proto-GCs had to be more 
massive than what is observed today (|Conrovll2011t IP'Ercole et al.ll2008l ). However, we show that multiple episodes 
of star formation can in fact take place during the lifetime of a cluster in particular between ^ 1 and ~ 10 Gyrs, 
thus suggesting a sizable enhancement in star formation. The overlap of this time range with the AGB phase further 
strengthens the case for stellar wind retention as a critical component in the formation of subsequent generations of 
stars in GCs. 

We thank Mark Krumholz, Charlie Conroy, David Pooley, Rebecca Bernstein, Chung-Pei Ma and Morgan Macleod for 
useful discussions. The software used in this work was in part developed by the DOE-supported ASCI/Alliance Center 
for Astrophysical Thermonuclear Flashes at the University of Chicago. Computations were performed on the Pleaides 
UCSC computer cluster. This work is supported by NSF: AST-0847563 (J.N. and E.R.), NASA: NNX08AL41G (J.N. 
and D.L.), and The David and Lucilc Packard Foundation (E.R.). 
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Fig. 1. — Average cluster mass loss rates and wind velocities as a function of time for Z = 1/10Zq. Green lines assume gas dynamics 
are dominated by the wind properties of the turn off stars, black lines show the population averaged values. Three representative times in 
the clusters age are denoted by the solid vertical lines. The current age of M15 is denoted by the vertical dashed line. 
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Fig. 3. — Hydrodynamic profiles for M15 with Mc = 4.5 X 10^ Mq and (t„ = Ukm/s fMcNamara et al.''2004; Gcrsscn ct al. 2002') at 
the three representative times in Figure [T] including the predicted profiles for its current age (dashed lines). Here, the population averaged 
prescription is used. 



0.0 



1.6 



01 
CD 
< 



105 j 



104 



103 



102 



IQi 



= 10^1 



10 



3.2 



4.8 



6.4 0.0 



1.5 



2.8 



4.3 



I 




20 



30 



Ov in l<m s" 



40 



50 5.0 



log(Mc) in Mg 



5.8 

CD 




0.16 



1.02 



2.47 



6.25 



15.85 



c 

— T 

o 

01 

5' 

® 
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Fig. 5. — Mass accumulation, Mace, as a function of cr for a fixed cluster core mass of Mc = 10^ Mq, calculated using the population 
averaged prescription. The blue shaded region delimits the boundaries for which our models collapse and trigger star formation. 
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Fig. 6. — Cycles of mass accumulation in a potential with _A/c = 10^ Mq, cr„ = 26km/s and ti = 2000 Myrs including star formation, and 
supernovae explosions. Mass accumulates till star formation is triggered, then 90% of this mass is assumed to form stars, removing the gas 
from the simulation (dashed lines). Supernovae persist for 100 Myrs at a rate of 2 X lO^^yr"-"^ further stripping mass from the system. 
For simplicity, we have assumed the mass injection properties to be steady during the simulation. 
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Fig. 7. — Gas properties for Ale = lO"^ Mq, Gv = 30km/s cluster at ti = 313 Myrs for different metallicities. Here, the population 
averaged values of the wind parameters are used. 
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